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Abstract 

We calculate the correction-to-scaling exponent lot that characterizes the 
approach to the scaling limit in multicomponent polymer solutions. A di- 
rect Monte Carlo determination of lot in a system of interacting self-avoiding 
walks gives lot = 0.415 ± 0.020. A field-theory analysis based on five- and 
six-loop perturbative series leads to lot = 0.41 ± 0.04. We also verify the 
renormalization-group predictions for the scaling behavior close to the ideal- 
mixing point. 

PACS: 61.25.Hq, 82.35.Lr, 05.10.Cc 



I. INTRODUCTION 



The behavior of dilute or semidilute solutions of long polymers have been investigated 
at length by using the renormalization group, 1-3 which has explained the scaling behavior 
observed in these systems and has provided quantitative predictions that become exact when 
the degree of polymerization becomes infinite. Most of the work has been devoted to binary 
systems, i.e. to solutions of one polymer species in a solvent. The method, however, can 
be extended to multicomponent polymer systems, i.e. to solutions of several chemically 
different polymers. The general theory has been worked out in detail in Refs. 4-6. In the 
good-solvent regime in which polymers are swollen, the scaling limit does not change. For 
instance, the radius of gyration R g increases as N", where 7 ' 8 v « 0.5876 and N is the length 
of the polymer. However, the presence of chemically different polymers gives rise to new 
scaling corrections in quantities that are related to the polymer-polymer interaction. In the 
dilute regime one may consider, for instance, the second virial coefficient B 2 between two 
polymers of different species. Its scaling behavior is 4-6 

B 2 = A(R g>1 R g , 2 )^ [1 + a(R gtl R g , 2 y^ + ■■■ + b{R gtl R g>2 )~^ + •••], (1.1) 

where R 9t \ and R 9t2 are the gyration radii of the two polymers, A, a, b are functions of 
Rg,i/ Rg,2, and ut, oj are correction-to-scaling exponents. Eq. (1.1) is valid for long polymers 
in the good-solvent regime; more precisely, for Aq, N 2 — > oo at fixed N\/N 2 (or, equivalently, 
at fixed Rg t i/R 9t2 ), where Aq and iV2 are the lengths of the two polymers. The exponent 
uj is the one that controls the scaling corrections in binary systems. The most accurate 
estimate of the exponent so far yields 7 A = uov = 0.517 ± 0.0071q;ooq, corresponding tow = 
0.880 ± 0.012l°;^ 7 . The exponent ujt is a new exponent that characterizes multicomponent 
systems. Perturbative calculations 9 ' 6 indicate that ujt is quite small, uot ~ 0.4. Thus, scaling 
corrections decrease very slowly in multicomponent systems and can be quite relevant for 
the values of Aq and N 2 that can be attained in practice. Therefore, the determination of 
the scaling behavior in multicomponent systems may require extrapolations in Aq and N 2 , 
which, in turn, require a precise knowledge of the scaling exponents. 
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In this paper we improve the previous determinations 9,6 of lot- First, we extend the 
perturbative three-loop calculations of Ref. 6. We analyze the five-loop expansion of ojt 
in powers of d = 4 — e, d being the space dimension, and the six-loop expansion of ujt 
in the fixed-dimension massive zero-momentum (MZM) scheme. Second, we compute ujt 
by numerical simulations. For this purpose we consider interacting self-avoiding walks and 
compute the second virial coefficient. A careful analysis of its scaling behavior provides us 
with a estimate of ujt- We also consider the ideal-mixing point where the interaction between 
the two different chemical species vanishes. A renormalization-group analysis of the behavior 
close to this point was presented in Ref. 6. An extensive Monte Carlo simulation allows us 
to verify the theoretical predictions. 

The paper is organized as follows. In Sec. II we present our renormalization-group cal- 
culations. In Sec. Ill we determine the correction-to-scaling exponent by means of a Monte 
Carlo simulation. In Sec. IV we discuss the ideal-mixing point where the effective interaction 
between the two chemically different species vanishes. Conclusions are presented in Sec. V. 

II. PERTURBATIVE DETERMINATION OF uj t 

The starting point of the calculation is the Landau- Ginzburg- Wilson Hamiltonian 6 

H= jd d x{\ [(d^) 2 + (d^ 2 f + n# + r 2 <g\ + 

^ o 0? + 2wo0?0! + ^]}, (2.1) 

where <f>i and 02 are n-component fields. As usual, the polymer theory is obtained in the limit 
n — > 0. In this specific case, the fixed-point structure of the theory is particularly simple 
and is explained in detail in Ref. 6. One finds that the (3 functions satisfy the following 
properties: (3 u (u,v,w) = ]3(u), (3 v (u,v,w) = ]3(v), (3 w (g,g,g) = 0(g), where (3(g) is the (3 
function in the vector 0(n = 0) y? 4 model and u,v,w are renormalized four-point couplings 
normalized so that u m Cuq, v ~ Cvq, w ~ Cwq at tree level. The relevant fixed point is the 
symmetric one u* = v * = w* = g*, where g* is the zero of (3(g). The exponent oot is given 
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by 



d(3 w 



(2.2) 



ijJT = 



dw 



The exponent lot can be computed directly in the 0(n — 0) ip 4 model. Indeed, as discussed in 
Ref. 6, ujt = —Va, where y 4 is the renormalization-group dimension of (fi^l in the symmetric 
theory with uq = Vq = wq. It corresponds to an 0(2n) vector theory and, for n — > 0, one 
is back to the 0(n = 0) (f 4 model. Using the results of Ref. 10 one can also show that, for 
n — ► 0, 0^02 is a spin-4 perturbation of the 0(2n) model and thus y& is the renormalization- 
group dimension of the cubic-symmetric perturbation Yla^t °f the 0(n = 0) ip 4 model. 
Thus, one can use the perturbative expansions reported in Refs. 11 and 12. The e expansion 
of ujt is 



At order e 3 it agrees with that given in Ref. 6. In the fixed- dimension MZM scheme we have 
at six loops 



where g is the four-point zero- momentum renormalized coupling normalized so that /3(g) = 



0.02. In order to obtain quantitative predictions, the perturbative series must be properly 
resummed. We use here the conformal-mapping method 13,19 that takes into account the 
large-order behavior of the perturbative series. 

From the standard e expansion we obtain in three dimensions (e = 1) uj t = 0.42 ± 0.04, 
while in the fixed- dimension MZM scheme we find u T = 0.37 ± 0.04. Using the pseudo-e 
expansion, Ref. 20 obtained cut = 0.380 ± 0.018 in the MZM scheme. Though compati- 
ble, the MZM result is lower than the e-expansion one. This phenomenon also occurs for 
other exponents and is probably related to the nonanalyticity of the renormalization-group 
functions at the fixed point. 14 ' 21 ' 15 ' 16 




(2.3) 




9 + Q 2 + 0(g 3 ), as used in, e.g., Ref. 13; the fixed point corresponds to 



14-18 



g* = 1.40 ± 
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A more precise estimate is obtained by considering ( = cot — co/2. The perturbative 
expansion of ( has smaller coefficients than that of uj t and thus ( can be determined more 
precisely. Its e expansion is 

C = — e 2 - 0.133936e 3 + 0.490572e 4 - 2.41405e 5 + 0(e 6 ) . (2.5) 
32 

The term proportional to e is missing, while the other coefficients are smaller by a factor of 
5-10 approximately. Similar cancellations occur in the MZM scheme: 

1 1 85 9 o 
C = "2 + 2 9 "«2 9 + 1368239 

-0.110394/ + 0.074425/ + 0.024718/ + 0(g 7 ). (2.6) 

Resumming the perturbative series, we obtain 

( = -0.006 ±0.009 (MZM), 

( = -0.008 ±0.012 (eexp). (2.7) 

We can combine these estimates with those for co. If we use the Monte Carlo result of Ref. 7 
reported in the introduction, we obtain 

co T = 0.433 ± 0.016lg;gg§. (2.8) 

If instead we use the field-theory estimates of cu reported in Ref. 17, we obtain 

UT = 0.399 ± 0.018 (MZM), 

cot = 0.407 ± 0.022 (eexp). (2.9) 
Collecting results, we estimate 

co T = 0.41 ± 0.04, (2.10) 
where the error should be quite conservative. 
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III. MONTE CARLO RESULTS 



In order to determine ujt numerically, we consider lattice self-avoiding walks (SAWs) 
with an attractive interaction — e between nonbonded nearest-neighbor pairs. If (3 = e/kT is 
the reduced inverse temperature, this model describes a polymer in a good solvent as long 
as (3 < fy, where /3$ ~ 0.269 corresponds to the collapse 9 transition. 22 ' 23 We consider two 
walks with different interaction energies ei and e 2 , i.e. with different (3i and (3 2 . We assume 
0i, @2 < Po, so that both walks are in the good-solvent regime. Then, we consider the second 
virial coefficient 

B 2 {Ni,N 2 -(3i,(3 2 ,(3i 2 )= l - J rf 3 r(l-e-^ 2 )) , r , (3.1) 

where the statistical average is over all pairs of SAWs such that the first one starts at the 
origin, has Ni steps, and corresponds to an inverse reduced temperature fii\ the second 
one starts at r, has N 2 steps, and corresponds to an inverse reduced temperature j3 2 . Here 
#(1,2) is the reduced interaction energy: if (1,2) = +oo if the two walks intersect each 
other; otherwise, H(l, 2) = —/3i 2 Af nn c, where Af nnc is the number of nearest-neighbor contacts 
between the two walks and 0i 2 = €i 2 /kT is the reduced inverse temperature. In order to 
generate the walks we use the pivot algorithm 24-28 with a Metropolis test, while the second 
virial coefficient is determined by using the hit-or-miss algorithm discussed in Ref. 29. We 
study the invariant ratio 

A (N N ■ ft ft ft \- ^(iVi,^;/?!,/^,/^) , v 

A 2 {l\i, 1\ 2 , fa, fj 2 , Pll) - rp (AT . R\T> (AT ■ R M3/2 ' 

where R g (N; f3) is the radius of gyration. As we have already discussed, in the limit Ni,N 2 — > 
oo at R g (N 2 ; (3 2 )/R g (Ni, (3\) fixed, A 2 obeys a scaling law of the form 6 

M N M ,p 2M ^f[^M), (3.3) 

where f(x) is universal. This applies for 0i, (3 2 < (3$ and, as we shall see, for (3i 2 sufficiently 
small. Note that all the dependence on the inverse temperatures is encoded in a function 
of a single variable. Moreover, f(x) is also the scaling function associated with a polymer 
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solution made of two different types of polymers that have the same chemical composition 
(hence (3 X = (3 2 = (3 12 ) but different lengths. In that case R g (N 2 ; f3 2 )/R g (N 1 ; fc) = (TVa/Aq) 1 ', 
so that Eq. (3.3) implies that A 2 (N 1 ,N 2 ;P,P,P) = g(N 1 /N 2 ), with g(x) = f(x v ) universal. 
The function f(x) satisfies the condition f(x) = f(l/x) and de Gennes' relation 30 f(x) ~ x p , 
p = 3/4-l/(2i/), iorx^O. 

We will be interested here in the corrections to Eq. (3.3). In the scaling limit we can 
write 

A 2 (N 1 ,N 2 ;(3 u( 3 2 ,(3 12 ) = f( P ) + ^ anm ltt. M f-(p) + ■ ■ ■ (3-4) 

n+m>l 

where p = R g (N 2 ; (3 2 )/R g (N l ] ft), x = [R g (N 2 ; p 2 )R g {N 1 ; A)] 1/2 , and we have neglected the 
contributions of the additional correction-to-scaling operators with renormalization-group 
dimensions —uj-i. They give rise to additional corrections proportional to x~ p , p = nux + 
muj + n i^i- Little is known about Ui, though we expect them to satisfy Ui > cut, > 
In the following we will assume that all such exponents satisfy uJi > 3ut ~ oj t + uj. The 
scaling functions g nm (p) are universal once a specific normalization has been chosen. Instead, 
the coefficients a nm depend on the model and, in particular, on the specific values of the 
parameters (3 2 , and P± 2 . 

We have simulated two SAWs with f3\ = 0.05 and (3 2 = 0.15, two values that are well 
within the good-solvent region. Then, we have computed A 2 for 100 < N\ — N 2 < 64000 
and several values of j3\ 2 in the range < fi 12 < 0.30. The results are plotted in Fig. 1. For 
f3 12 < 0.25, as iV = Ni = N 2 increases, the estimates of A 2 tend to become independent of 
P12 although the convergence is very slow. The behavior changes for (3 12 > 0.25 and indeed 
the data indicate that A 2 = for iV — > oo for some fi 12 = p i2tC slightly larger than 0.25. This 
value corresponds to the case in which the short- distance repulsion is exactly balanced by the 
solvent-induced attraction proportional to fi\ 2 . In field-theoretical terms, this means that 
the renormalization-group flow is no longer attracted by the symmetric fixed point discussed 
in Sec. II, but rather by the unstable fixed point with w* = 0. Thus, at f3\ 2 = (3i 2:C there 
is effectively no interaction between the chemically different polymers. For p i2 > (3i2, c , A 2 

7 



becomes negative signalling demixing. The behavior of A 2 does not change significantly if f3\ 
and f3 2 are varied. In Fig. 2 we report results for shorter walks for different pairs of Pi and 
P 2) and also for walks with Ni = AN 2 . Note that Pi 2>c depends very little on the parameters 
Pi and p 2 . 

In order to determine lot, we have considered the data with p 12 = 0, 0.05, 0.10, 0.15 that 
are sufficiently far from the critical value Pi 2jC . We performed a fit that is linear in uj t 

HMN; P12 = 0) - A 2 {N- Pi 2 )\ = ci(Pi 2 ) -cu T lnx+ + (3.5) 

where Pi 2 = 0.05, 0.10, 0.15. The exponents uo a and uo b should take into account the additional 
scaling corrections. Since 2u T ~ uj, we should have oo a xs uo T ph uo — uo t . Using the field-theory 
estimate of u T and to « 0.88±0.03 (Ref. 7), it should be safe to take 31 to a = 0.43±0.06. As for 
uo b we should have uo b m 2uo T a;. Moreover, there is also the possibility that there exists an 
additional correction exponent U\ not very much different from 3u>t, which would contribute 
a correction with exponent ui — lot- For this reason we have taken uj b = 0.85 ± 0.20. The 
error should be large enough to include all possibilities. The results are reported in Table I. 
The systematic error reported there gives the variation of the estimate as u a and uj b vary 
within the reported errors. The results with iV m i n = 250 and 500 are compatible within 
errors and thus we can take the estimate that corresponds to N niin = 500 as our final result. 
To be conservative, however, the error bar takes also into account the possibility that the 
observed small trend is a real one. If the neglected corrections are of order x~ 3uJT we expect 
the results to depend on N mhl as N~?™ T pa iV~j n 7 . This implies that the estimate of u T can 
decrease at most by 0.005 when 7V min is further decreased. This leads to the result 

cut = 0.415 ± 0.020. (3.6) 

As a check we perform a nonlinear fit of the form (fit 2) 

AmM =A t + ^ + ^ + ^ (3.7) 

since, x ~ A^, we can also fit the data to (fit 3) 



A 2 (N;P 12 ) = A* 2 + 



ai{Pu) a 2 (/3 12 ) a 3 (/3 12 ) 



(3.8) 



pjA T j\T2A T jy3A T ' 



where At = lotv. The results are reported in Table I. They agree with those obtained 
before and allow us to estimate the universal constant A* 2 [A 2 = f(p) for p 1.24; f(p) is 
defined in Eq. (3.3)]: A* 2 = 5.495 ± 0.020. 

IV. IDEAL-MIXING POINT 

In this section we consider the behavior close to the ideal-mixing point (IMP) /3 12)C 
where the effective interaction between the two chemically different species vanishes. The 
renormalization-group analysis is presented in Ref. 6. The behavior is controlled by an unsta- 
ble fixed point characterized by an unstable direction with renormalization-group dimension 
yi = 2/u — 3 and by a stable direction with exponent — u, where uj m 0.88 is the usual 
correction-to-scaling exponent. These results imply that close to the IMP a renormalization- 
group invariant quantity 1Z (for instance, the invariant ratio A 2 introduced above) scales 
as 



In this section we wish to verify this scaling behavior for the second virial coefficient. For 
this purpose we have made simulations for six different pairs of f3\ and (5 2 with N\ = N 2 
and 0.25 < /3 12 < 0.28. Since N x = N 2 = N and R g {N; (3) w a{(3)N u , we can rewrite the 
previous equation as 



K((3 12 ) = fn[(Pi2 - Pi 2 , c )(R 9 ,iR 9 , 2 ) yi/2 ,R 9 ,i/R 9 , 2 ] 

+ (» I ^ /2 9n[{Pi2- (3i2, c ){R 9 +R 9 , 2 T'\ R 9 ,i/R 9 ,2}. 




(4.1) 



lliN-^fc,^) = f n (b,p) + -^gn&p), 
b = (/3i2 - Pn,c)Nt 



p = R g (N;P 1 )/R g (N;P 2 ), 



(4.2) 



where (3\ 2tC also depends on (3\ and (3 2 , and 
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= vyi = 2 - 3u = 0.2372 ± 0.0003. (4.3) 

The critical value (5 12 ,c can be characterized by requiring A 2 (N — > oo; f3 2 , Pi2,c) — 0. We 
can also define a finite-iV IMP as the value of /? i2 where A 2 (N; (3i, ft 2 , P12) vanishes (this is 
the analogous of the Boyle point in 9 solutions): we define Pf 2iC (N) such that 

A 2 (N;P 1 ,p 2 ,pf 2jC (N)) = 0. (4.4) 

Inserting in Eq. (4.2) we obtain for N — > 00 the behavior 

Pll(N) = (3 12 , c + j^. (4.5) 

Finally, we can replace fi 12:C with P^ C (N) in Eq. (4.2) obtaining the equivalent form 

U(N; p u p 2 , p 12 ) = f n [(p 12 - 0*JN+, p] + ^^[(A2 - tf&N+, p], (4.6) 

where g n {b, p) vanishes at the IMP 6 = 0. Eq. (4.6) is more suitable for a numerical check 
close to the IMP than Eq. (4.2), since scaling corrections vanish at the IMP and are therefore 
small close to it. In the following we verify numerically predictions (4.5) and (4.6). 

In Fig. 3 we show f3l 2 c (N) vs N'^'^ for three different pairs of p 1 and (3 2 . The data 
show a quite good linear behavior: only the two points corresponding to iV = 2000 and 4000 
are in some cases off the linear fit, probably because our sampling is not yet adequate for 
these large values of N. These results allow us to obtain /3i 2)C : 



Pl2,c = 


0.2574 ± 0.0006 


A 


= 0.05, fa 


= 0.10, 


Pl2,c = 


0.2588 ±0.0007 


A 


= 0.05, fa 


= 0.15, 


Pl2,c = 


0.2609 ± 0.0003 


Pi 


= 0.05, (3 2 


= 0.20, 


Pl2,c — 


0.2596 ± 0.0004 




= 0.10, p 2 


= 0.15, 


Pl2,c = 


0.2609 ± 0.0009 


Pi 


= 0.10, (3 2 


= 0.20, 


Pl2,c = 


0.2626 ± 0.0009 


Pi 


= 0.15, fa 


= 0.20. 



Note that the dependence on (3\ and /3 2 is tiny. 
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Then, we verify Eq. (4.6). In Fig. 4 we report A<i vs the scaling variable {fi^i — Puc)^^ 
for different Pi, p 2 . All points with 100 < N < 1000 fall on top of each other confirming the 
scaling behavior predicted by the renormalization group. The scaling functions depend on 
Pi and p 2 through the ratio p. Moreover, one should also take into account that the scaling 
variable is only defined up to an arbitrary prefactor. We now show that the dependence on p 
is tiny for the range of values of p we have considered, since all data with different values of 
N, Pi, and p 2 fall on a single curve once one takes as scaling variable R(Pi, p 2 )(Pi 2 — Pl 2c )N ct ', 
where R(Pi, p 2 ) is a properly chosen constant that depends on Pi and P2. This is evident 
in Fig. 5, where we report data with different Pi and p 2 . We only consider N = 500 for 
clarity, since we have already verified that data with different values of N show the predicted 
scaling. The scaling is very good, indicating that the p dependence is negligible. If we choose 
i?(0.05, 0.10) = 1 the scaling curve can be parametrized as A 2 = —19.773a; — 46.457x 2 , with 
x = R(Pi, p 2 )(Pi2 — Pnc)^- Finally, note that R(Pi,p 2 ) does not depend very much on Pi 
and p 2 . For instance, i?(0.15, 0.20)/i2(0.05, 0.10) « 0.88. 

V. CONCLUSIONS 

In this paper we have considered the corrections to scaling expected in multicomponent 
polymer solutions. A high-precision Monte Carlo simulation with very long walks, up to 
N = 64000, gives uo T = 0.415 ± 0.020. This estimate is consistent with the estimate ujt = 
0.41 ±0.04 obtained by using the perturbative renormalization group. Previous perturbative 
renormalization-group calculations gave oo T ~ 0.37 (Ref. 9) and uo T w 0.40 (Ref. 6): they 
substantially agree with our estimate. On the other hand, the numerical result 6 ujt ~ 0.35 
obtained by exploiting the relation between ojt and the growth exponent for four-arm star 
polymers seems slightly too small. 

One should note that ujt is quite small and thus convergence may be quite slow. For 
instance, in the case we have considered numerically, A* 2 = 5.495 ± 0.020 for iV — > 00. 
On the other hand, for /3 12 = (resp. P12 = 0.15) we find A 2 = 5.790 ± 0.005 (resp. 
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A 2 = 5.232 ±0.005) for N — N± — N 2 — 64000. Even if the walks are very long, there is still 
a 5% discrepancy. For N = 1000, differences are larger, approximately of 15%. 

In order to obtain a better qualitative understanding of the corrections we have also per- 
formed additional simulations. The results are reported in Fig. 2. The qualitative behavior 
is very similar in all cases and almost independent of fa and (3 2 . In particular, corrections 
appear to vanish in all cases for 0.05 < f3 12 < 0.10 and to increase strongly for @ 12 > 0.20. 
Morever, (3\ 2fi is always close to f3$, and shows a tiny dependence on (3\ and (3 2 . 

The discussion presented here addressed the behavior of multicomponent solutions, but 
it should be noted that the results are also relevant for copolymers in which chemically 
different polymers are linked together. 32 Also in this case scaling corrections with exponent 
ujt are present. 

Finally, we have also considered the behavior close to the ideal- mixing point, where 
there is no effective interaction between the two chemically different polymer species. Our 
numerical results are in very good agreement with the renormalization-group predictions of 
Ref. 6. 
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3045 (2005), since the two polymers are ideal. In this case scaling corrections are controlled 
by a new exponent uj$, whose e expansion is u>s = e — \e 2 + 1.15154e 3 — 4.20439e 4 + 
16.9338e 5 + 0(e 6 ). A Pade-Borel resummation gives uo s = 0.77±0.02. A three-loop analysis 
gives u s ~ 0.82 (Ref. 6). 
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TABLES 



type 


N ■ 




A* 


x 2 


DOF 


fit 1 


100 


0.424 ± 0.002 ± 0.004 




15.2 


20 


fit 1 


250 


0.418 ± 0.004 ± 0.006 




12.4 


17 


fit 1 


500 


0.415 ± 0.007 ± 0.008 




11.8 


14 


fit 2 


100 


0.426 ± 0.015 


5.498 ±0.011 


32.2 


26 


fit 2 


250 


0.421 ± 0.030 


5.495 ±0.018 


32.1 


22 


fit 3 


100 


0.426 ± 0.015 


5.499 ±0.011 


32.1 


26 


fit 3 


250 


0.423 ± 0.029 


5.496 ± 0.017 


32.0 


22 



TABLE I. Results of the fits. When two errors are quoted, the first one is the statistical error, 
while the second one is the systematic error. DOF is the the number of degrees of freedom of the fit. 
Fit 3 provides an estimate of At = vojt- We compute lot = At/v by using v = 0.5876 ± 0.0001. 
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FIG. 1. Invariant ratio A 2 for ft = 0.05, ft = 0.15 vs ft 2 for several N = N\ = N 2 . The two 
figures differ only by the vertical and horizontal scales. Lines connecting points with different ft 2 
are only intended to guide the eye. 
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FIG. 3. Finite- N ideal-mixing point /3ff c vs l/N w+< ^ for three different pairs of /?i and The 
line corresponds to the fit /3ff c = /?i 2)C + a/TV^+f 
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FIG. 4. Invariant ratio A2 vs (/?i2 
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FIG. 5. Invariant ratio Ai vs R{fin — Pi2, c )N^ for N = 500. We report results corresponding 
to six different pairs of Pi and 
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